CMB polarization from secondary vector and tensor modes 



Silvia Mollerach 
Centro Atomico Bariloche, Av. Bustillo 9500 
8400 Bariloche, Argentin^ 

Diego Harari 

Departamento de Fisica, FCEyN, Universidad de Buenos Aires 
Ciudad Universitaria - Pab. 1, 14-28, Buenos Aires, Argentin^\ 

Sabino Matarrese 
Dipartimento di Fisica 'G. Galilei', Universita di Padova 
INFN - Sezione di Padova 
via Marzolo 8, F35131 Padova, Italffl 
(Dated: February 2, 2008) 

We consider a novel contribution to the polarization of the Cosmic Microwave Background in- 
duced by vector and tensor modes generated by the non-linear evolution of primordial scalar per- 
turbations. Our calculation is based on relativistic second-order perturbation theory and allows to 
estimate the effects of these secondary modes on the polarization angular power-spectra. We show 
that a non- vanishing B-mode polarization unavoidably arises from pure scalar initial perturbations, 
thus limiting our ability to detect the signature of primordial gravitational waves generated during 
inflation. This secondary effect dominates over that of primordial tensors for an inflationary tensor- 
to-scalar ratio r < 10 -6 . The magnitude of the effect is smaller than the contamination produced 
by the conversion of polarization of type E into type B, by weak gravitational lensing. However the 
lensing signal can be cleaned, making the secondary modes discussed here the actual background 
limiting the detection of small amplitude primordial gravitational waves. 

PACS numbers: 98.70.Vc,98.80.Es,98.80.Cq 



I. INTRODUCTION 



The generation of a stochastic background of primordial gravitational waves is a fundamental prediction of inflation- 
ary models for the early Universe. Its amplitude is determined by the energy scale of inflation, which can widely vary 
between different inflationary models. The detection of this gravitational radiation would provide a crucial test for the 
validity of the whole scenario. Gravitational wave detectors, however, are quite unlikely to have enough sensitivity to 
detect such a primordial signal, owing both to its smallness and to its extremely low characteristic frequencies. The 
existence of ultra-low-frequency gravitational radiation, however, can be indirectly probed thanks to the temperature 
anisotropy and polarization it induces on the Cosmic Microwave Background (CMB) radiation. In particular, the curl 
component, called B-mode, of the CMB polarization provides a unique opportunity to disentangle the effect of tensor 
(gravitational- wave) from scalar perturbations, as this is only excited by either tensor or vector modes 0,0- From this 
point of view, future satellite missions, such as Planck, which will have enough sensitivity to either detect or constrain 
the B-mode CMB polarization predicted by the simplest inflationary models, might represent the first 'space-based 
gravitational- wave detector' The main complication in this context comes from the effect of gravitational lensing 
on the CMB by the matter distribution, which implies the transformation of E-mode into B-mode polarization Q: 
such a non-linear effect might actually obscure the signal due to primordial tensor modes. It has been pointed out 
that the inflationary gravitational-wave background can only be detected by CMB polarization measurements if the 
tensor to scalar ratio r > 10~ 4 , which corresponds to an energy scale of inflation larger than 3 x 10 15 GeV [HISSIS- 
Quite recently, however, a better technique to clean polarization maps from the lensing effect has been proposed, 
which would allow tensor-to-scalar ratios as low as 10 -6 , or even smaller, to be probed 0,^3- Other secondary 
contributions to the B-t ype polarization arising during the reionization stage, though of much smaller amplitude, 
have been considered in fill . 
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The case of vector modes is even more interesting, as they cannot be produced during inflation and are in general 
extremely difficult to generate at early epochs, with the exception of models which predict the existence of primordial 
tangled magnetic fields 0, 0] . 

This paper considers a new source of B-mode polarization, coming from secondary vector and tensor modes. The 
contribution to temperature anisotropy arising from these modes has already been analyzed in Refs. |l4l Il5l|. The 
evolution of cosmological perturbations away from the linear regime is in fact characterized by mode-mixing, which 
not only implies that different Fourier modes influence each other, but also that primordial density perturbations act 
as a source for curl vector perturbations and gravitational waves |43( [lfj. Let us emphasize that these secondary 
vector and tensor modes always exist and that their amplitude has a one-to-one relation with the level of density 
perturbations, which is severely constrained by both CMB anisotropy measurements and Large-Scale Structure obser- 
vations. Therefore, their properties are largely inflation model-independent, contrary to primary tensor modes whose 
amplitude is not only model-dependent, but is well-known to be suppressed in some cases, like e.g. in the so-called 
curvaton model for the generation of curvature perturbations ^3 • 

The plan of the paper is as follows. In Section II we introduce the second-order vector and tensor modes which are 
produced by the non-linear evolution of primordial scalar perturbations. In Section III we obtain the contribution of 
these secondary modes to the polarization angular power-spectra, while in Section IV we compare these contributions 
to those from primordial gravitational waves and gravitational lensing. Section V contains our main conclusions. 



II. SECOND-ORDER VECTOR AND TENSOR MODES 



The perturbed line-element around a spatially flat FRW background takes a particularly simple form in the so-called 
Poisson gauge [T^! which to linear order reduces to the Newtonian gauge. Adopting the conformal time dr\ = dt/a, 
one can write 

ds 2 = a 2 (r]) {- (1 + 2*) drf - 2V i dr ] dx i + [(1 - 2$) % + 2H tj ] dx l dx j } , (1) 

with 5ij the Kronecker symbol. The lapse perturbation is the sum of a first-order or primary term (indicated by a 
label P) and a second-order one (indicated by a label S): ^ = <3>p + ^s- 

The shift perturbation Vi in this gauge is a pure vector, i.e. d 1 Vi = 0, which only arises as a second (or higher)-order 
contribution (i.e. Vi = Vsi). The spatial metric perturbations contain a scalar mode, which includes both a linear 
and a second-order term, $ = $ P + $ s , and a tensor (i.e. transverse and traceless) mode {d l Hij = H\ i = 0), 
which also contains first and second-order contributions, namely Hij = Hpij +Hsij- Here and in what follows spatial 
indices are raised by the Kronecker symbol 8j . 

Hereafter we assume that the Universe is spatially flat and filled with a cosmological constant A and a pressureless 
fluid - made of Cold Dark Matter (CDM) plus baryons - whose stress energy-tensor reads T^ v = p u^u u (u^u^ = — 1). 
This is a reasonable aproximation for the purposes of this paper, although the inclusion of the radiation in the 
evolution equations would give more accurate results at the smaller scales considered, that entered the horizon during 
the radiation dominated era. 

In this paper we will concentrate on vector and tensor perturbations. In the standard inflationary scenario for the 
origin of perturbations linear vector modes are not present, but they are generated after horizon crossing as non- linear 
combinations of primordial scalar modes. 

Let us report here the second-order calculation which leads to the generation of these secondary vector and tensor 
modes. For the A = case, this problem was originally solved in Ref. starting from the results of second-order 
calculations in the synchronous gauge and transforming them to the Poisson gauge, by means of second-order 

gauge transformations |2l|. Here we will solve the problem directly in the Poisson gauge and we will account for a 
non- vanishing A term (see also p^ - ). 

The background Friedmann equations read 3H 2 = a 2 (8irGp + A) and p — —3H.p, where dots indicate differentiation 
with respect to r/, Tt = a/ a and p is the mean mass-density. Also useful is the relation 2"H = —Ti 2 + a 2 A. 

We then perturb the mass-density and fluid four- velocity as p = p{\ + 6) and u M = (<5q + v^)/a, with 5 = <5p + <5g 
and = Vp + v§ . 

Let us briefly report the results of linear perturbation theory in this gauge (for a detailed analysis, see e.g. [23jp. 
The non-diagonal components of i-j Einstein's equations imply ^Pp = $p = cp, for the scalar part, and 

H Fj + 2HH Pj - V 2 iJ Pj - , (2) 

for the tensor part, while its trace gives the evolution equation for the linear scalar potential ip, namely 



(p + 'iUip + a 2 kip = 



(3) 
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Selecting only the growing-mode solution, we can write ip(x,rj) = ipo{x.)g(r]), where ipo is the peculiar gravitational 
potential linearly extrapolated to the present time (?7o) and g = D + /a is the so-called growth-suppression factor, 
where D + (rf) is the linear growing-mode of density fluctuations in the Newtonian limit and a the scale-factor. In the 
A = case g = 1. An excellent approximation for g as a function of redshift z is given in Refs. 0,113 



5«O m [nfl 7 - fi A + (1 + fim/2) (1 + flA/70)] , (4) 

1/2 

with a m = Q Qm (l + z) 3 /E 2 (z), n A = n aA /E 2 (z), E(z) = [rj 0m (l + z) 3 + n 0A ] and Qorn, Qoa = 1 - O 0m , the 
present-day density parameters of non-relativistic matter and cosmological constant, respectively. We will normalize 
the growth-suppression factor so that g(z = 0) = 1. 

The energy and momentum constraints provide the density and velocity fluctuations in terms of ip, namely 

' [VV - 3H {ip + Hip)] , 



AwGa 2 p 

vpi = — A i di{<p + Htp) . (5) 

Perturbations of the matter stress-energy tensor up to second-order in the Poisson gauge have been calculated in 
Ref. |2f|, for a general perfect fluid. Specializing to the pressureless case, one has 

T° = -p(l + 6 P + 6 s +v P ) , 

T* = -p[v P +vi + (ip + 5 P )v P ] , 

T ,. = P~[vpi + vsi-Vi + (-3ip + 6p)vpi] , 

T) = p v P v Pj , (6) 

where v P = vpv P j. Note that the second-order velocity vsi is the sum of an irrotational component uS , which is the 

gradient of a scalar, and a rotational vector vjt , which has zero divergence, d l v^ = 0. 

For the purpose of obtaining secondary vector and tensor modes, we can start by writing the second-order 
momentum-conservation equation |44j . — 0, which gives 

(y$ —Vij+Ti (v£) - Vi) = -d&s - tyvpi - \d, (vp + ip 2 ) . (7) 

For pure growing-mode initial conditions ip oc ip and vpi oc diip, which makes the RHS of this equation the gradient of 

a scalar quantity; thus, the vector part only contains a decaying solution (vg^ — Vi) oc a^ 1 , and we can safely assume 

v {1) -V- 
v Si — **• 

To proceed one needs the second-order perturbations of the Einstein tensor, S^'G^: these can be found for any 
gauge in Appendix A of Ref. [27j ■ 

The second-order 'momentum constraint' S^G l = iirG5^T l gives 

d l (Hf s + <f>s) + '^V 2 V i + ipd i ip + ±<pd i ip = -AnGa 2 p (<p + 5 P ) v P + . (8) 

The pure vector part of this equation can be isolated by first taking its divergence to solve for the combination 
Tt^S + and then replacing it in the original equation. We obtain 

V 2 V 2 V % = 16nGa 2 p d J (v Pj d t 5 P - vp l d j S P ) . (9) 

We can further simplify this equation and write: 

W 2 V Z = ~F{z) (d^oVVo - &&<pod m + 2fl,-e ) , (10) 

with 

= 2g 2 (z)E(z)f(Q m ) ^ 



4 



where [24],|25| /(0 TO ) = dhiD + /dh\a w il m (2:) 4 / 7 , i/o is the Hubble constant and 



v 2 e = -~ ((vVo) 2 - d i d kVo d i d k VQ ) 



(ii) 



For A = 0, the above expression for Vi reduces to Eq. (6.8) of Ref. [Tg, noting that F(z) = r\ in that case. 
Finally, we can solve for the second-order tensor modes by looking at the traceless part of i-j Einstein's equations, 
8^G l } - \5^G k k 5) = SnGiS^T) - \8 { ^T k k 5)), namely 



iv 2 ($ s - *s) + | (V^) 2 + -tpV 2 ? 



d l d 3 (<J> S - s ) + 28 i <pd j tp + Aip&djtp 



~ (&Vj + djV*) + H (d% + djV*) + (S* aj + ^nH^ - V 2 I%) = 8nGa 2 p (v P v Pj - ±v 2 P sfj 



(12) 



To deal with this equation, we first apply the operator $di, which allows to solve for the combination $s — ^S; 
and then replace it in the original equation, together with the expression for the vector mode Vi . After a lengthy but 
straightforward calculation we obtain 

(13) 



V 2 V 2 (W sj + 2HH l S] - V 2 ff^) = V 2 d k d t R t k 5) + 2V 2 (V 2 ff j - d l d k R k J - d k djR\) + d l d j d k d i R i k , 

where we introduced the traceless tensor 

= &<pdH<p - \ (Vipf Si + AnGa 2 p (v P v Pk - 

2 / 2E 2 (z)f 2 (n m ) \ ( t l m , 2 xi 

= 9 l 1 + 30ow(1 + , ) 3 1 (^0^0-3^0) S k 



R> 



(14) 



Equation (|13|) can be solved by Green's method, as the corresponding homogeneous equation is the one for linear 
tensor modes, whose analytical solutions are known in the limiting cases O m — * 1 or J1a — * 1. In the A — ► limit one 
recovers the result of Ref. 0] , namely 



Hsi 3 (x, rj) 



(2tt) 3 



d 3 fce ,k.x40 5y(k) 



1 _ jijky) 
3 krj 



where jl are spherical Bessel functions of order I and iSiy(k) = J d 3 ye y Sij(y), with 

Sij = V 2 6 5ij + didjOo + 2 (<9i<9,c/? V 2 (/3o - did k (p d k djip ) 



(15) 



(16) 



For the purpose of the present analysis, the main contribution comes from the early evolution of Hgij, when the 
cosmological constant term can be neglected. One can then adopt the Einstein-de Sitter result above, up to an overall 
correction factor g 2 ^ = g 2 (z — > oo). 

It proves convenient to express the previous quantities in Fourier space. Let us define 



1 



(2tt) s 



d 6 k 



(k, 77)Qj +1) (k, x) + V ( -^ (k, r,)Ql-V (k, x) 



where (2^ 



( 3j ±1) ( k ' x ) = ^Tf^ 1 ±ie 2 )jexp(ik-x) , 



with ei, e 2 and 63 = k forming an orthonormal basis and 

V^(K V ) = ^^ J i^(e lTi e 2 ).k' (fc 2 -2k.k')^o(k'Vo(k-k'). 

The peculiar gravitational potential tpo is a Gaussian random field with Fourier space correlation 

(^o(k)^o(k')) = (2^) 3 P)(k + k')/yfc) , 



(17) 



(18) 



(19) 



(20) 



5 




where P v (k) is the gravitational potential power-spectrum P v (k) — P 0ip k^ 3 (k/ko) nB ~ 1 T 2 (k), k is some pivot wave- 
number and T s is the usual matter transfer function [30l |. which is unity on large scales and drops off like k~ 2 on small 
scales because of the stagnation effect of matter perturbations on sub-horizon scales during radiation dominance. For 
later convenience we can relate P v to the dimensionlcss power-spectrum of the comoving curvature perturbation (see 
e.g. HH) K, namely A 2 K (k) = A 2 , (/coX^Ao)™^ 1 , where A 2 n {k ) = (25/9)(P 0v /27r 2 ) 0. 

The vectors V^ 1 ^ are two independent, chi-square distributed random fields each one with power-spectrum 



Pv(k, v ) 



1 F(z) 2 



9ir 2 k 4 

which, after numerical integration yields 

36tt 2 



/l pOO 
dcos6>(l-cos 2 6>) J dk'k' 4 (k 2 - 2kk' cos6>) 2 P v (k!)P^ k n + k 2 — 2k' k 



Pv{k) 



C v (n s )A 4 7l (ko)(kF(z)) 2 k- 



W v {k/K) , 



cos (9) , (21) 



(22) 



(25) 2 " v-z-vcv-m- ^ 

where fc* = 17o m /i 2 Mpc -1 , with ft, the Hubble constant iJo in units of 100 km s _1 Mpc -1 ; the coefficient CV(ng) 
ranges from 0.062 to 0.29 for n s between 1 and 0.8, respectively. The function Wy(x) (which is also weakly dependent 
on n s ) is unity for i«l and drops to zero at x « 1; its form is plotted in Figure 1 (for n s = 1). A good fit, useful 
for next section calculations, is given by Wy{x) = (1 + 5x + 3a; 2 ) -5 / 2 . 
For tensor modes we similarly define 



i 



(2tt) 3 



H {+ V (k, V )Q\+ 2) (k, x) + (k, ry)Qj- 2) (k, x 



(23) 



where |2£ 




Q^ ±2) (k,x) = -\/-(ei ± ie 2 )i(ei ± ie 2 )j cxp(zk • x) 



The primary tensor modes can be expressed as 



H p ±2 \k lV )=x P " 



(±2) n ,?>ji{kri) 



j (k) 



(24) 



(25) 



for the Einstcin-de Sitter case. 



spectrum P XP (k) 



where x^p 2 ^ 



( — 2) 

and Xp are two independent Gaussian fields each one with power- 
\k/k Q ) nt T 2 {k). The tensor transfer function T t (k) reads T t (k) = [1 + l.Uy + 2.5y 2 } 1 / 2 
[32j, where y = k/k cq and /c cq = 0.073flo m h 2 Mpc -1 is the horizon scale at matter radiation equality. In the A ^ 
case the evolution of primary tensor modes is more complicated and will be dealt with using CMBFAST |3^| . 



Pq X p k~ 
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We can similarly relate P XP (k) to the dimensionless tensor power-spectrum (see e.g. A 2 (fc) = ^i(ko)(k/ko) nt , 



where A 2 (fc ) = (Poxp/^k 2 )- Defining the 'tensor-to-scalar' ratio as in Ref. [3l||, we find 

^ 24 ^T- ( 26 ) 

In single-field slow-roll inflation models the standard consistency relation n t — — r/8 also applies [34j . 
The secondary tensor modes are instead characterized by 

ff r 2) (k^)«(i-^)^xr 2) (k), (27) 

where we have introduced the factor g 2 ^ that makes this a good approximation in the A ^ case for z > few as 
discussed above. The secondary tensors in Eq. I|27l) are given by 

x ^\ k ) = -^— J |^_[ k '.(e lT e 2 )] 2 ^(k'Vo(k-k'), (28) 
and are chi-square distributed independent random fields with power-spectra 

rj r -i p 1 pOQ 

p xs(. k ) = zf^2tf J d cos (9 sin 4 (9 y dk'k' 6 P v {k')P v {^k' 2 + k 2 - 2k' k cos 0) , (29) 
which can be numerically integrated to give 

Pxs = ^C xs AUko)k- 3 fj-j (^J W xs (k/K) , (30) 

where the coefficient C xs (ns) ranges from 0.062 to 0.29 for n s between 1 and 0.8, respectively. The function W xs (x) 
(which weakly depends on n s ) is unity for i < 1 and drops to zero at x sa 1; its form is also plotted in Figure 1 (for 
n s = 1). It can be well fitted by W x (x) = (1 + 7x + 5x 2 ) -3 . 



III. POLARIZATION ANGULAR POWER-SPECTRA 



Polarization of the CMB arises from Thomson scattering of anisotropic radiation by free electrons. The generation 
of temperature and polarization anisotropics in the CMB from gravitational perturbations has been studied in detail 
by several authors [J, 0, l29l l35j | . A simple and powerful formalism is the total angular momentum method |29l l35| , 
which includes the effect of scalar, vector and tensor modes on an equal footing. This is the most convenient approach 
for our calculations, and we are going to extensively use the results of Ref. 29]. The temperature and polarization 
fluctuations are expanded in normal modes that take into account the dependence on both the angular direction of 
photon propagation n and the spatial position x, s G™(x, n) 

3 1 m T 2 (31) 

(Q±iu)( V , x ,n) = f E (^ (m) ±^ m) ) ±2 Gr, 

J \ ' l ro=-2 

with spin s = describing the temperature fluctuation and s = ±2 describing the polarization tensor and m = 
0, ±1,±2 denoting scalar, vector and tensor perturbations, respectively. and are the angular moments of 

the electric and magnetic polarization components and 

s Gr(x,n) = (-0y^L*nn)]ex P (ik.x). (32) 

The Boltzmann equation describing the time evolution of the radiation distribution under gravitation and scattering 
processes can be written as a set of evolution equations for the angular moments of the temperature, Oj (for £ > to), 
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and both polarization types, E^ and Bf 1 ' (for £ > 2 and m > 0) 



(m) 



(m) 







K l Q M _ 0^1+1 qW 



(2^-1) (2^ + 3) 



e+i 



-E 



2m 



(21-1) l - 1 t{i + i) 



B 



(m) _ 2 K £+l g(m) 



(2^ + 3)" 
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(m) 



2m 



(2£- l)"'- 1 
where the coupling coefficients are 

The fluctuation sources are given by 



-E. 



(m) _ 2^1+1 ^g(m) 



£(^ + 1) 4 (2* + 3) 



e+i 



rB [ r\ 



{P-m 2 ){P~s 2 ) 



( 2 



(33) 
(34) 
(35) 

(36) 



S (2) = f p(2) _ £ ^ 



(37) 



with 



1 r 



p(m) = J_ 

10 



(38) 



The modes with m = — \m\ satisfy the same equations with £^~' m '^ = — B^" and all the other quantities unchanged. 

These equations can be formally integrated, leading to simple expressions in terms of an integral along the line-of- 
sight 29| . The temperature fluctuations are given by 



21 + 



f- = / d V 'e-^S ( ^(i)jf m \H^-v')) 



(39) 



where jf are given in ref. |29j . For the polarization, we have 



21+ 1 



(40) 



where the radial functions read 



,(±i) 



(x) = -!)(£ + 2) 



-M*)+i^) + 2#+4^ 
a;^ a; 



j' e (x) + 2 



(41) 



(42) 



Here the differential optical depth f = n e /7T& sets the collision rate in conformal time, with n e the free electron density 
and (Jt the Thomson cross section and t(t)o, 77) = f t(r]')dr)' the optical depth between r\ and the present time. The 
combination fe~ T is the visibility function and expresses the probability that a photon last scattered between dr\ of r\ 
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and hence is sharply peaked at the last scattering epoch. In early reionization models, a second peak is also present 
at more recent times. 

Scalar modes do not contribute to B-polarization, thus — 0. We are interested in the contribution to the 
angular power-spectrum for the E and B modes arising from vector (m = 1) and tensor (m = 2) perturbations, 

, (B) _ 2 rdk l3 |E«'">(,„,m 2 



7T J k ^ (2£+ l) 2 ' 



A. Effects of decoupling 

Before recombination, photons, electrons and baryons behave as a single tightly coupled fluid, so we can find 
approximate expressions for the polarization source 

pirn) 

using a perturbative expansion in inverse powers of the 
differential optical depth f (fc/f <§C 1) (3(|. For the vector modes we need to consider the Euler equation for the 
velocity perturbation of baryons to close the system 

4 ) =V-%% ) -V) + UqU- v U ) . (44) 
a it 

To leading order in the tight coupling approximation we obtain 

BP=v%\ (45) 

and from the baryon Euler equation 

v ( s ] ~ V. (46) 

To this perturbative order the quadrupole vanishes. However it does not vanish to the next order, and thus also 
^ 0. An appropriate expression for our calculation can be obtained by combining the equations for and 
E^' to obtain an evolution equation for and replacing in it the leading order solution we get 

p(D + L+pW _ Jl^y = o. (47) 
10 10V3 



This can be integrated to give 



P«(M) = ^=J\'V{k,rf)e X p{~T{ V ',r,)). (48) 



Then, the B-mode polarization multipoles can be estimated by 



2^+1 10 



driTe- T (3^(k(vo-v)) £ d V 'V(k,ri')exp (-^T(r)',r,)\ , (49) 



An analogous expression holds for Eg replacing by ep . Because of the sharpness of the visibility function 

around the time of decoupling r)r>, V(k, rj') and /3g(k(r)o — r))) (or ) can be evaluated at r)r> and taken out of the 
integration. The remaining integral can be evaluated analytically approximating the visibility function by a Gaussian 
of width A??n pi 

B \ e { l°f = -^O.SlAr^Mz^^o - no)). (50) 
Then, the resulting angular power-spectrum of B modes from vector perturbations results 



Cf m = ^l 2 ^ 2 D I dk k A [/# 3 (fc(r? - Vd)) Pv(k,r} D ), (51 ! 



2 



9 



where we have added the contribution coming from the m — ±1 modes. For ci^^ we obtain an analogous expression 



with ef'\k(rio — Vd)) replacing ft^ (k(r)Q — fyo)). Replacing Eq. l{2"2"|) we obtain 



t{i+l)C\ 



B)(l) 



2tt 



= £{i+l) 16^A^(fc )(fe 



2(n 3 -l) 



Arw) ( F{z D ) 



'la 



Vok* 



(52) 



The integral can be conveniently evaluated using the WKB approximation for the Bessel functions appearing in /3 
and 4 , leading to 



(i) 



<(4 1} ) 2 (*)) ~ + 2) 



1 



x^\J x 2 — I 2 



(53) 
(54) 



for x > £, and for x < I. The results for the vector contribution to the B and E angular power-spectrum are 
shown in Figure 2, where we have used the following values for the parameters, for the concordance ACDM model 
(fi A = 0.7, Q m = 0.3, tt b h 2 = 0.022, h = 0.71): Ar) D /r) D = Az D /2z D = 0.09, rfy/r) D = 30, K = 0.135 Mpc^ 1 and 
A^(fco) = 2.3 x 10~ 9 at fc = 0.002 Mpc" 1 . 

For the tensor modes, the tight coupling approximation to the leading order gives P^ — /4 = —H^/Sf. To 
the next order we obtain the following evolution equation for P^ 



o rr(2) 
p(2) + ± fP (2) + ^_ = 

10 10 

that can be integrated 

P {2) (k, V) =-Jof d? ^ (2) ( fc ' ^) expC-^r^', r,)). 
Proceeding as for the vector case, we obtain for the B modes 

and the same result changing f3 t by e\ for the E modes. 
Using Eq. (Sfll we can write 



(55) 



(56) 



(57) 



Hf(k,r lD )=-Zg\ 



32{kr) D ) (2 ) 



>)D 



x^(fc). 



Then, for the angular power-spectrum, we obtain 



t^ + i)cf ){2) 

2^ 



£(^+1)^0.51 



Aqc 
VD 



<& I dk k 2 W\k( m - VD ))} P xs (k)jl(kr) D ), 



(58) 



(59) 



for the B modes (and the same expression with /3» 2 ' — ► ep for the E modes), with P Xs (k) given by Eq. H30(l . To 
perform the remaining integration it is convenient to use the WKB approximated expressions 



<($ 2) ) a (s)> 
{{eff{x)) 



l / Vx 2 - £ 2 



x 3 Vx 2 — i 2 

2 



(60) 



1 - 



l + £(£ + l)/2 



1 



Vx 2 - f 2 



x\J x 2 — i 2 



for x > £, and for x < £. The results are shown in Figure 2. The dominant contribution comes from the tensor 
modes for low multipoles, I < 20, while the vector contribution dominates for larger I. 
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B. Effects of reionization 



The recent detection by WMAP of excess power on large angular scales in the correlation between temperature and 
E-polarization of the CMB is a strong evidence that the Universe was reionized at relatively early times [37], u9 ■ The 
best fit to the Thomson scattering optical depth due to reionization is given by r ri = 0.17, with still large uncertainties. 
In a spatially-flat ACDM model with £Ia — 0.7, reionization would have happened at a redshift z r ; = 17, if it took 
place in a single step. 

Early reionization provides another opportunity for Thomson scattering to modify the polarization properties of the 
CMB. Its effects can be approximately decomposed in two parts [3!| . On the one hand, rescattering damps already 
present anisotropics and polarization by a factor e~ T ". On the other hand, rescattering of the existing quadrupolc 
anisotropy significantly enhances the polarization signal at low multipoles. 

Here we estimate the effects of reionization on the polarization signal induced by secondary vector and tensor 
modes. We shall assume, for simplicity, that reionization took place in a single step. We can then approximate in 
Eq. H40f> that the visibility function fe~ r is sharply peaked at Tfo and thus evaluate /3i m \k(r)o — -q)) (and ) at 
X] = ?7ii and take them out of the integral. The coupling between photons and electrons after reionization is, however, 
not very tight, and thus we can not proceed with the same approximations as in the previous section to compute the 
polarization sources p( m \r)). Indeed, f = —OTn e a with n e — 0.88nto(l + z) 3 X(z) and X(z) the ionization fraction. 
Thus r = 0.0019-ffo(l + z) 2 X(z)(ilbh 2 / 0.022) (0.71/ h), and for the cosmological parameters under consideration, 
k/f w 0.5fc77o at z = z r \, and the tight coupling approximation is not applicable at any relevant wavelength. 

The polarization produced after reionization is essentially due to the scattering of existing quadrupole anisotropies, 
generated by scattering at earlier times or by gravitational effects during propagation. Since the quadrupole polariza- 
tion Ei generated during decoupling is typically smaller than the temperature quadrupole O2, we can approximate the 
polarization source as p( m ) « O^/IO, and use for the formal solution of Eq. H39[). dropping in the fluctuation 
sources S of Eq. (|37[1 the terms proportional to f against V (in the case of vectors) or against H (for the tensors), 
since the coupling is not tight. 

The contribution of secondary vectors to B-polarization after reionization is then given by 

- - -— ^dTTT-e-^fcfoo -V)) dv'eM-riv'^MKr,')^ 1 ^(fcfa- 77')), (61) 



2^+1 2 .,, (n 

where 32 \x) — y/3j2(x)/x. For a single step reionization scenario, we can approximate 

We have defined an effective "reionization width" A?/ r ; = rj T i f Arjie~ T ', which for the TO = 0.3, spatially-flat ACDM 
model with single-step reionization and optical depth r r j = 0.17 is Af^-i = 0.16?7 r i. This can be calculated taking into 
account that, for z < z r ;, 

Notice that r r ; = 0.17 for z r ; = 17 if fi m = 0.3. The reionization contribution to the power-spectrum of B-polarization 
by vector modes is finally approximately given by 

"' lK " ;1 -W + l) (I) 2 A^fco) (fe)^"^ Cy(n s ) (^) 2 (Z^y ( % - (64) 



2tt 



where y = xrj T i/(r]o — rjn). An analogous expression holds of course for C\ E)(1) with ef ] in place of (3^. 

The contribution of secondary tensor modes to the polarization induced after reionization can be estimated with 
similar approximations: 

^U+'l^ = ~1T f°^re- T 4 2) (k( Vo -v)) J\v^M-r{i,v))H (2 \k^)^ 2 \k{ 1 -i)) 1 (65) 
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FIG. 2: Angular power-spectrum of B (left panel) and E (right panel) polarization by secondary vector (solid lines) and tensor 
modes (dashed lines), in a ACDM model without reionization (thin lines) and with an optical depth to reionization th = 0.17 
(thick lines). 



where j^ 2 \x) = 3 j2{x)/x 2 . We approximate 

B \ e { l°f = -^§A Vli (l + S)H^(k, Vl ,)(3 { e 2 \k( m TH))D(krH) • (66) 

We have defined an additional coefficient S such that SAr] T i = d?yf e _T J v At]' e~ T{ - r) >' q >. It accounts for the fact that 
in this case (opposite to the vector case) the time dependence of the source after reionization is not negligible. For 
the model parameters used throughout this work to estimate the effects, 6 ~ 0.4. The damping factor D(x) = 4/x if 

(22) 

x < 4 (and D = 1 otherwise) arises because j\ (x) can be approximated as 1/5 only for small x; larger values lead 
to oscillations that make the integral decay as 1/fc. Finally, the reionization contribution to the power-spectrum of 
B-polarization by secondary tensor modes is approximately given by 



l(£+ l)cf )(2) 27 /Ar^ 2 



- ^+^^2 [ — ) (! + Wo / dfc fe 2 $Hk(m - fM)) P Xs (k)m^)D''(kVn) ■ (67) 
2tt 2bir z \ rjii J J L J 

The numerical result for the angular power-spectrum of E and B polarization by secondary vector and tensor modes 
produced after reionization, in the concordance ACDM model, is displayed in Fig. 2, added to the polarization that 
was produced during decoupling, damped by a factor e~ 2TVi . 



n 2 



IV. COMPARISON WITH PRIMARY TENSOR MODES AND COSMIC SHEAR 



The signal from the secondary modes computed in the last section must be compared with the signal from the pri- 
mordial gravitational-wave background generated during inflation. This depends on the model of inflation considered 
and its amplitude is parameterized by the 'tensor-to-scalar' ratio r (see Eq. (HHJ)- Its angular power-spectrum can 
be accurately computed by available codes like CMBFAST against which we have tested the accuracy of our 
semi-analytic approximations in the previous sections. The polarization induced by primordial gravitational waves can 
be calculated as in Eq. ifCTl) replacing the power-spectrum of secondary modes P xs (k) by that of the primary modes 
P XP (k). The approximate expressions, both with or without reionization, are in good agreement (with more than 20% 
accuracy) with the numerical results obtained using CMBFAST. Notice that the polarization induced by primordial 
gra vitational waves may be slightly smaller than it was calculated here due to damping by neutrino free-streaming 
[40| . that we did not take into account. 
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FIG. 3: Angular power-spectrum of secondary total (vector plus tensor) B-polarization (thick solid line) compared with that 
induced by primordial gravitational waves (dashed line) with a tensor to scalar ratio r = 10~ 6 , in a flat ACDM model (fl m = 0.3) 
with Tri = 0.17. The thin solid line is the signal due to gravitational lensing cleaned by a factor 40. 



In Figure 3 we plot the numerical result (calculated with CMBFAST) for the B-polarization induced by primordial 
gravitational waves in a model with tensor to scalar ratio r = 10 -6 , along with the total secondary B-modes (vector 
plus tensor) derived in this work. It is clear that the effects of primary and secondary modes are comparable. Primary 
modes have somewhat larger effects at the lowest (I < 10) multipoles, and also for I between 30 and 100, but secondary 
modes overcome them in intermediate and larger values of I. The power-spectrum induced by primordial gravitational 
waves scales as r. Thus, secondary vector and tensor fluctuations limit the ability to detect primordial gravitational 
waves through measurements of B-polarization if r < 10~ 6 . 

The extraction of the B-polarization signature of primordial gravitational waves is seriously limited by weak gravi- 
tational lensing effects, that convert the dominant E-type polarization into B-modes pj, which swamp the signature 
of primordial tensor fluctuations unless they have a ratio to scalar fluctuations r of the order or larger than 10~ 4 . 
Fortunately, the weak lensing effect can be accounted for, through reconstruction of the gravitational lensing potential 
by the correlations of B-polarization between large and small angular scales, which primordial gravitational waves do 
not produce. "Cleaning" of the gravitational lensing signature may achieve factors of 40 in the power-spectrum, or 
even larger [Tol |. In Figure 3 we display the predicted gravitational lensing signal in B-polarization calculated with 
CMBFAST for the cosmological model under consideration, reduced by a factor 40. Clearly, an improvement in the 
cleaning algorithm would convert the secondary vector and tensor modes derived here in the barrier to the detection 
of primordial gravitational waves through B-polarization of the CMB if r < 10~ 6 . 



V. CONCLUSIONS 



The study of the magnetic-mode polarization of the CMB will become a fundamental and possibly unique tool to 
search for the stochastic gravitational- wave background, whose detection would represent a clear signature of a period 
of inflation in the early Universe. In contrast with temperature anisotropies and E-type polarization of the CMB, 
scalar perturbations do not give rise to B polarization in a direct way, so that measurements of the B-mode could 
be used to probe rather small gravitational-wave background amplitudes. Because of this reason, much observational 
effort is taking place for its detection, and future dedicated missions, such as NASA's Beyond Einstein Inflation Probe 
|4l| or ground-based experiments, like BICEP an d PolarBeaR, are being planned. 

The main background for the detection of the B-mode is represented by the gravitational lensing conversion of a 
fraction of the dominant E-type into B-type polarization. However, it has been recently shown 10] that the lensing 
signal can largely be cleaned, thus allowing to probe inflationary models with tensor-to-scalar ratios r < 10 -6 . As the 
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largest signal compared to the lensing one comes from low multipoles, where the reionization contribution is dominant, 
this limit depends on the Thomson scattering optical depth to reionization, that still has a large uncertainty. 

At these sensitivity levels, there are other secondary effects that can give rise to sizable contributions to the B-type 
polarization. We estimated here the contribution coming from secondary vector and tensor modes, which originate 
by the mildly non- linear evolution of primordial density perturbations. The amplitude and harmonic content of this 
contribution is completely fixed, once the primordial power-spectrum of the density perturbations is known. For 
a concordance ACDM model, and adopting the high reionization redshift implied by the WMAP data [oTI l3S| , we 
found that the effect of secondary vectors and tensors becomes comparable to that of primary gravitational waves for 
r < 10 -6 , thus making it the actual background for the detection of primordial gravitational waves through B-mode 
polarization. 
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